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Abstract 



We present a short review of theoretical results (mainly for ex- 
perimentalists) published in many different papers. The formulae are 
^ . presented for the different integrated cross sections, the number of 

interacting nucleons, multiplicities of secondaries, the dispersions of 
multiplicity distributions. Two possible tests for the search of Quark- 
tCT* . Gluon Plasma formation are discussed. CERN SPS data for pro- 

duction of secondaris in central Pb-|-Pb collisions are compared with 
Quark-Gluon String Model predictions. 
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1 Introduction 

At present time, different Monte Carlo codes are used, as a rule, for the 
analyses of existing heavy ion experiments and for planning the future ones. 
However, sometimes it is not clear enough, what physics was used in some 
Monte Carlo code, which processes were accounted for, and so on. In the 
present paper we give a review of the simplest theoretical predictions which 
can be used for different estimates, not involving complicated calculations. 

Formally we consider the high energy nucleus-nucleus collision as a super- 
position of the independent nucleon-nucleon interactions. However, really 
the main part of the discussed results is based practically only on geometry, 
and does not depend on the model of interaction. 

2 Cross sections 

Let us consider the collision of nucleus A with nucleus B, assuming that 
they are not very light. The simplest expression for the total inelastic 0"^^, 
or secondary production cross section a^^^ comes from geometry: 

a'^^^aZ' = 7riRA + RBr. (1) 

The difference between a'^J^^ and o"^*^ , which corresponds to all processes 
of disintegration or exitation of one or both nuclei without secondary pro- 
duction, is negligibly small at high energies and heavy incident nuclei. In 
the case of heavy ion collisions the total elastic cross section is only slightly 
smaller than a'^J^^ [0], af^ ^ <^^ab-> so 

<B = <B + ^AB = MRa + Rb? ■ (2) 

More accurate is Bradt-Peters expression |0] 

a'XB = T:Rl{A^I'' + B^'''-c)\ (3) 

which accounts for the possibility of very peripheral collisions without in- 
elastic interactions. The values of parameters Rq and c are slightly different 
for light and heavy nuclei 0]. The comparison of this expression with the 
experimental data was presented in 0. 

More detailed results for different cross sections of heavy ion interactions 
can be obtained with the help of multiple scattering (Glauber) theory, see 
[|] for details. This theory can be used at all energies starting from several 
hundred MeV per nucleon. At high energies (say, more than 10 GeV per 
nucleon) in principle the inelastic screening [^ ^ should be accounted for. 



however in the case of heavy ion cross sections these effects are numerically 
very small. 

The amplitude oi A — B elastic scattering with momentum transfer q can 
be written in the frame where i?-nucleus is a fixed target as 

FUq) = '^ J d'he^'^'ll - SABib)] , (4) 

where k is the incident momentum of one nucleon in A-nucleus, b - impact 
parameter, and 

SAB{b) = {A\{B\Y[[l-TNN{b + u,-s,)] \B)\A) , (5) 



jeB 



where Ui and Si are the transverse coordinates of nucleons. Contrary to the 
case of hadron-nucleus interaction, we can not integrate Eq. (5) analitically 
even after standard assumption that nuclear density is a product of one- 
nucleon densities. 

To make the problem manageable, one can pick up a certain fraction 
only of various contributions in the expansion of the product in Eq. (5), 
distinguished by large combinatorial factors. The leading graphs correspond 
to the so-called optical approximation [0, which sums up the contributions 
with scattering not more than once for each nucleon, i.e. takes into account 
the products of amphtudes Ftvat in Eq.(5) with all indices i, j being different. 
It corresponds to the summation of diagrams shown in Fig. 1 and it contain 
the largest combinatorial factor A{A — l)...{A — n + l)B{B — l)...{B — n + l) in 
the case of n-fold interaction. In all another series the combinatorial factors 
are significantly smaller, however these series have different signs, so due to 
the large cancelation some of them can give significant contribution to the 
final expression. To avoid crowding of lines we have only shown in Fig. 1 
the nucleon participants from the nucleus A (upper dots) and B (lower dots) 
with the links standing for interacting amplitudes, and we do not plot the 
nucleon-sp ect at or s . 

In the language of Eq. (5) the optical approximation can be reproduced, 
making the averaging (74|...|y4) and {B\...\B) inside the product Hij iii Eq.(5) 
instead of averaging the product as a whole: 

S7Bib) = Ei^mil - TMb + u. - s,)] \B)\A) , (7) 

Using the standard assumptions of the multiple scattering theory one 
obtains 

SjBib) = [l - jT^Ab)]^ ^ exp[-T,,,{b)] (8) 



with 

T,pt{b) = a/2J<fhTA{h - h)TB{h) (9) 

and 

Ta(6) = A / dzp{h = Vr2 - z^, z) (10) 

J —oo 

where p{r) is the one-particle nuclear density, J p{r)d^r = 1. 

In the Reggeon language the optical approximation corresponds to the 
accounting only one pole (nuclear ground state) in the both M^ and M^ 
complex planes instead of all intermediate states, see [Q. 

However, numerical calculations ^ (see also ^ for the case of light nuclei 
collision) demonstrate that the optical approximation is not accurate enough 
even for the total cross sections. The difference with the data amounts ~ 
10-15% for 0"^^ and is much larger for differential cross sections. 

The more explicit rigid target [jlO|, [Tl|, |12|, |l3l approximation can be re- 
produced, making the averaging only {B\...\B) inside the product in Eq. (5) 

SAUb) = {A\{Yl\{B\[l-TMb + u^-s,)] \B)}\A) = [T,.,.(&)]^ (H) 






Tr.gXb) = J J d\TA{h - b)exp 



^Tsibi] 



(12) 



This approach corresponds to the diagrams of Fig. 2, where each nucleon 
from the nucleus A can interact several times, but all interacting nucleons 
from B are still different. Due to the obvious asymmetry in contributions 
of the two nuclei such approach can be justified theoretically in the limit 
A/B <^ 1, say for C — Pb, ot S — U collisions, however sometimes it can be 
used for the case of heavy ion collisions with equal atomic weights, see next 
Sect. 

Further corrections to the elastic amplitude have been considered in Refs. 



14, 15, 16, 17, 18|, however these results are rather complicated for practical 



use. 



Another possibility is the direct calculation of Eq. (5) using Monte Carlo 
simulation [0, |l^ . 

Let us present now the list of several integrated cross sections. The total 
cross section oi A — B collision is given by the optical theorem 



tot 



2 J d%[l ~ SAB{b)] , (13) 

where the function SAsib) can be taken, say, from Eq. (8) or Eq. (11). The 
integrated cross section of elastic A — B scattering is 

a% = ld''b[l-SAB{b)f. (14) 



The sum of cross sections of all diffractive processes without production of 
secondaries, i.e. the elastic and quasielastic scattering can be calculated 
with the help of close approximation J2a' \^'){^'\ = 1 similarly to the case 



of hadron- nucleus scattering |20 



'^AB = J dM^ - '^SABib) + [lAB{b)f} , (15) 

where in the optical approximation 

iTBib) = 1 - \RoAb) , (16) 

Roptib) = a'^'' J (fbiTAibi - b)TBibi) (17) 

and 0"*"*^' is the A^A^ inelastic cross section. In the rigid target approximation 

lABib) = ^jd%,TA{b, - 6)exp[-a*"^'TB(6i)] (18) 

The correspondent expression for the secondary production cross section 
reads 

<B = <B - <yAB = I dMl - [IaB^} ■ (19) 

It is also possible [^ to obtain the integrated stripping cross sections of 
collisions among light nuclei. This can be done by considering the different 
intermediate states P^ , p3| , |^ of all Glauber diagrams for elastic AB am- 



plitude. The resulting expressions have the simple form of a combination 
of several total inelastic cross sections. For example, the cross section for 
stripping one nucleon from nucleus A in A-B collision, cr^^, is written as 

^S = ^«l - ^(7-i)b) , (20) 

where the density distribution of the nucleus with atomic weight A-1 should 



be slightly "corrected" [21 



3 Distributions on the number of interacting 
nucleons 

Let us consider the events with secondary hadron production in nuclei A and 
B minimum bias collisions. In this case the average number of inelastically 
interacting nucleons of a nucleus A is equal pi] to 



A prod 

< Na >....= ^ . (21) 

(^ab 



If both nuclei, A and B are heavy enough, the production cross sections 
of nucleon-nucleus and nucleus-nucleus collisions can be written as 

aTB=^Rl. aTB' = <RA + RB)\ (22) 

and accounting for Eq. (1) we obtain that in the case of equal nuclei, A = B, 
for minimum bias (m.b.) events which are averaged over impact parameter 

< Na >m.b= A/A =< Na >c /4 . (23) 

So in the case of minimum bias events the average number of interacting 
nucleons should be four times smaller than in the case of central collisions 



25], where (neglecting by small corrections which will be discussed below) 
< Na >c~ a. 

The average number of inelastic A^A^ interactions in minimum bias AB 
collisions with any secondary production is equal 0, |T^ to 

< i^AB >m.b= ABa'S'^^/aTB'' ■ (24) 

Let us note that the optical approximation cannot be used for the 
calculations of the averaged numbers of interacting nucleons and the distri- 
bution on this number. The rigid target approximation gives here reasonable 
results. In particular, for the distribution over the number of inelastically 
interacting nucleons A^^ of A nucleus in AB interaction we have |]T^ 



'''<^^' = l^ iA-N.m^ I ''^"l^-f")!'-"* II - ^-("'1"-^ • (25) 



(^AB 



{A-Na)\Na\ 



The theoretical predictions for the distributions on the electrical charge 
of non-interacting nucleons based on Eq. (25) were compared with the data 
in 



1^ and the agreement was quite reasonable, see Fig. 3 taken from [|n 
All the presented results are related to the case when only the events 
with secondary production are registrated. In the case when the events with 
nuclear disintegration without secondary production are also registrated, one 
should change a^f, by a^Jv and a^As^ by a'Xe- 

Eq. (25) is written for minimum bias events. In the case of events for some 
interval of impact parameter h values, the integration in Eq. (25) should be 
fulfilled by this interval, hmtn < b < hmax- In particular, in the case of central 
collisions the integration should be performed with the condition h < bo, 
and bo -C Ra- In the case of the collisions of equal heavy ions the value of 
< Na > decreases with increase of the impact parameter. As even at small 
6 7^ some regions of colliding ions are not overlapping, only very small 
fraction of events, not larger than 0.5-0.7% of all minimum bias sample, can 
be considered as the central interactions. 

6 



In the case of asymmetrical ion collisions, say S — U, all nucleons of light 
nucleus at small impact parameters go throw the regions of relatively high 
nuclear matter density of heavy nucleus, so practically all these nucleons 
interact inelastically. For the case oi S — U interactions at ^/SNN = 20 GeV 
all events with 6 < 2 ^ 3 f m can be considered as the central ones |^ . 

The predictions for the distributions on the number of inelastically interac- 
ting nucleons at different impact parameters can be found in [p5l . 

4 Multiplicities of secondaries, inclusive den- 
sities and tests for QGP formation 

Let us now look at the inclusive density of the produced secondaries, 

Pab{x) = E^^ . (26) 

In the case of asymptotically high energies we predict^ for xj? ^ |]^ 

Pab{x)/pnn{x) =< vab > (27) 

and for secondary multiplicities 

< nAB > / < ^NN > = < l^AB > (28) 

However at the realistic energies both right-hand side ratios (27) and 
(28) are significantly smaller than the predicted values due to the energy 
conservation corrections. For example, if one nucleon of the projectile nucleus 
interact with several nucleons of the target nucleus, the total energy of all 
secondaries should be equal to the initial energy, so the effective energy of 
every NN interaction is smaller than the initial energy. These effects are 
more important in the both A and B fragmentation regions. 

The corrections to Eqs. (27) and (28) decrease with the growth of the 
initial energy. Some predictions for the multiplicity of secondaries with ac- 
counting the division of energy between several NN interactions can be found 

The relations between average multiplicities and inclusive densities of sec- 
ondaries produced in heavy ion collisions can be used for search of the Quark- 
Gluon Plasma (QGP) formation. One of them is the relation between average 

^We do not account here the additional shadowing due to percolation [g6[ or multiple 
pomeron interactions p^. 



multiplicity and inclusive density in the central region. In the case of A^A^ 
collisions for any z = n^j^]^ / < ti^nn >) the ratio | 

R{z) = pnn{x)/ < Pnn{x) > U~o (29) 



is a non-trivial function of z p9|. The reason is again in energy conservation. 
Let us compare two events, first with ?t,_ =< n_ > and second with n_ = 2 < 
n_ >. Evidently, if the shapes of the inclusive spectra in these events will be 
the same, the energy of all secondaries in the second event will be two times 
larger than in the first one. So the shapes should be different, in the second 
event we will find smaller number of fast secondaries and larger number of 



slow ones. The experimental data [Q as well as the model calculations p9 
show that the ratio R{z) = Pnn{x)/ < Pnn{x) > \x^o is about 3.5 at z = 2.5. 

However, in the case of heavy ion collisions with independent nucleon- 
nucleon interactions we can not expect any difference of Eq. (29) from the 
linear function, because every independent nucleon-nucleon collision gives 
the same contribution to p(x) as well as to multiplicity 0. The possible viola- 
tion of such linear dependence should be considered as an evidence for some 
collective interaction (possibly, QGP formation). The analysis of the data 
at energy about 4 GeV per nucleon [^ shows good agreement with linear 
dependence, and the similar analysis of the data at higher energies seems 
to be very interesting. Such analysis can be provided in some kinematical 
domain. 

The second test for QGP formation is connected with analysis of the mul- 
tiplicities of different secondaries produced in the central and minimum bias 



A — B collisions ||2^. In the case of independent nucleon-nucleon collisions 



the multiplicity of secondaries produced in the central region should be pro- 
portional to the number of interacting nucleons of projectile nucleus, that 
is four times larger than in minimum bias events, see Eq. (23). It should 
depends also on the average number, < z/^vb >, of inelastic interactions of 



every projectile nucleon with the target nucleus. As was shown in |22[, the 
average number of interactions in the case of central nucleon-nucleus colli- 
sions, <v>ci is approximately 1.5 times larger than in the case of minimum 
bias nucleon-nucleus collisions, < i'>m.b.- After accounting several another 



^The considering ratio is rather similar to KNO one, however instead of cross sections 
we use inclusive densities. 

■^The difference in heavy ion and NN interactions here is connected both with experi- 
mental conditions and with confinement. In heavy ion collisions we usually registrate only 
the interacting nucleons and neglect the non-interacting ones. Due to that the energy of 
secondaries can change from event to event. In the case of NN interaction the spectator 
quarks and gluons will be involved into the interaction due to confinement forces, and the 
energy of all secondaries is exactly equal to the initial energy. 



corrections [^ one obtain 

<n>,~4.5 <n>^,b. ■ (30) 

In the conventional approach considered here, we obtain the prediction 
of Eq. (30) for any sort of secondaries including pions, kaons, J/i^, Drell- 
Yan pairs, direct photons, etc. Let us imagine that the QGP formation is 
possible (or dominate) only at comparatively small impact parameters, i.e. 
in the central interactions. In this case Eq. (30) can be strongly violated, 
say, for direct photons and, possibly, for light mass Drell-Yan pairs, due to 
the additional contribution to their multiplicity in the central events via the 
thermal mechanism. At the same time, Eq. (30) can be valid, say, for high 
mass Drell-Yan pairs and for pions, if the most part of them is produced 
at the late stage of the process, after decay of the plasma state. So the 
experimental confirmation of Eq. (30), say, for pions and its violation for the 
particles which can be emitted from the plasma state should be considered as 
a signal for QGP formation. Of course, the effects of final state interactions, 
etc. should be accounted for in such test. 

5 Multiplicity distributions and dispersions 



It was shown in Ref. |Q that the main contribution to the dispersion of 
multiplicity distribution in the case of heavy ion collisions comes from the 
dispersion in the number of nucleon-nucleon interactions. The last number 
is in strong correlation with the value of impact parameter. 

For the normalized dispersion D/ <n>, where D^ =<n^> — <n>^ we 
have [53 



<n>2 



where 
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(31) 



(32) 



is the average number of nucleon-nucleon interactions in nucleus-nucleus col- 
lision, n and d are the average multiplicity and the dispersion in one nucleon- 
nucleon collision. 

In the case of heavy ion collisions < pab >~ 10^ — 10^, so, as a rule, the 
second term in the right hand side of Eq. (31) becomes neghgible [Q, and the 



first term, which is the relative dispersion in the number of nucleon-nucleon 
interactions, dominates. In the case of minimum bias A-B interaction the last 
dispersion is comparatively large due to large dispersion in the distribution 
on A^^ {vab)- So in the case of some trigger (say, J /ip production) without 



fixing of tlie impact parameter, tlie multiplicity of all another secondaries 
can change significantly in comparison with its average value. 

In the case of some narrow region of impact parameters we have the oppo- 



site situation. Dispersion in the distribution on A^^ {^ab) is very small [^ 
especially in the case of central collisions. The dispersion in the number of 
inelastic interactions of one projectile nucleon with the target nucleus, i'nb, 
should be the same or slightly smaller in comparison with the minimum bias 
case. So the dispersion in the multiplicity of the produced secondaries can 
not be large. 

If fluctuations in the effective number of interactions are negligible, the 
second term in the right hand side of Eq. (31) should dominate and the 
dispersion of the produced secondaries will be suppressed by the factor < 
i^AB > in comparison with the case of A^A^ collisions. 

In the case of small fluctuations in u^b one can see that any trigger can 
not change significantly the average multiplicity of secondaries in the central 
heavy ion collisions, even if this trigger strongly influents on the multiplicity 
in the nucleon-nucleon interaction. 

Some numerical example of very narrow secondary multiplicity distribu- 
tion in central heavy ion collisions can be found in ||2^ 



Now let us consider some rare event C in high energy heavy ion collisions. 
It can be Drell-Yan heavy mass pair, J/\l/, T, W or Z production. The nor- 
mahzed multiplicity distribution of the produced secondaries in the events 
associated to such rare events, Pc{n), should be not the same as the distribu- 
tion in the mimimum bias events, Pm.b.i''^)-, with Y^ Pm.b.{n) = J2 Pc{n) = 1. 
The reason is quite clear. The rare events can occur independently in every 
A^A^ collision, so in the case of heavy ion interactions the probability of rare 
event is proportional to uab, Eq- (32). So ||3^ 



Pc{n) = -^^P{n) (33) 

< ^AB > 

The experimental data at high energies are in agreement with the predic- 
tion of Eq. (33). An example of such agreement is presented in Fig. 4, taken 
from 12^. The similar relation can be written for the distributions on the 



transverse energy, E-r |35, 36 1, etc., and they are in agreement with existing 



high energy data 0, ^ . 



6 Spectra of secondaries 

As we discussed, the energy conservation effects violate asymptotical ratios 
(27) and (28). The possible way to account for these effects was suggested 
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in pHl |3^. The inclusive spectrum of secondaries produced in nucleus A 



nucleus B collision can bew written as 



1^ dy - mR\(v)<N.>-^ ^.^ ' + (34) 

where y is the rapidity of secondary h in cm. frame and the functions 
^Asiy) account for the energy conservation effects. Here we connect the 
inclusive spectra of secondaries in the heavy ion collisions with the spec- 
tra in the nucleon-nucleus interactions. To calculate the last ones as well 
as the functions -R^b(2/) ^^ ^^e one of the variants of the dual topological 
unitarization approach - the quark-gluon string model (QGSM) |^. This 
model has been used successfully to describe the spectra of secondaries pro- 



duced on both nucleon [^, ^, ^] and nuclear |^, ^ targets. In QGSM 
the high-energy interaction is determined by an exchange of one or several 
pomerons. One pomeron is coupled to a valence quark and diquark, and 
the others are coupled to sea quark-antiquark pairs. All elastic and inelastic 
processes are the result of a cut of different number of pomerons or a cut 
between pomerons |^ . The inclusive spectra of secondaries are determined 
by convolutions of the quark and diquark momentum distributions and the 
functions describing their fragmentation into hadrons. All these functions are 
determined by their Regge asymptotics. When one nucleon interacts with 
several, the average number of pomerons is increased, which tends to soften 
the momentum distributions of all quarks and diquarks and, consequently, 
the secondary hadron spectra. 



The functions RAsiy) can be taken in the form 



< jy >NA fi-yA) 

where f^{y, < z/ >) is the contribution from the interaction with u nucleons 
to the spectrum of secondary particle h in NA collision. 

The results of secondary spectra calculation in the QGSM depend on sev- 
eral parameters. They can be fixed from pp data. The values of parameters 
were changed slightly in comparison with |3^, because new experimental 



data have appeared. The results for the x^^-spectra of secondary 7r+, vr , K^ 
and K~ are presented in Fig. 5. The agreement with the data [H^, E3| is 



quite reasonable. The description of secondary proton and antiproton spectra 



[agreement with the data is of the same order) one can find in [^ 



The rapidity distributions of secondary vr"*", vr , K^, K , p and p measured 



in Pb+Pb central collisions at 158 GeV/c per nucleon |^ are compared with 
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our calculations using Eq. (34) in Fig. 6 .In the case of pions the agreement 
is reasonable. The differences between the calculated curves and the data 
are larger than 10%. In the case of secondary kaons we see that there is a 
problems. 

The spectrum of K~ is in good agreement with the model, whereas for K^ 
disagreement is about 20%. In the case of change of the model parameters 
one will obtain the disagreement in Figs. 5c and 5d. Even more serious is 
the K^ / K~ ratios. Experimentally this ratio at y ~ is about 1.5, however 
in the model we can not obtain it more than 1.2. It means that the sea 
contribution (which leads to the same yields of K^ and K~) is too large in 
the model. It is possible, that the account for percolation/string fusion effect 
p^ , |T7| , ^ as well as final state interactions |^ can solve this problem. 



In the cases of p and p the agreement with the data is good enough. 

7 Conclusions 

We presented some simplest and model independent formulae obtained with 
assumption that high energy heavy ion collision can be consider as the super- 
position of independent nucleon-nucleon collisions. So any serious disagree- 
ment of these predictions with the data can be considered as the signal for 
some collective interaction, say QGP formation. 

Contrary to the case of elementary particle physics the numerical results 
for heavy ion collisions should be obtained correspondent^ to the experi- 
mental conditions. For example, if only events with secondary production 
are registrated one should use Eqs. (21) and (24), whereas if events with nu- 
clear disintegration and without secondary production are also registrated, 
the cross sections a^'^ , a^^ and a^^fj should be replaced by (J^abi '^nb ^^^ 
(J^i^%, respectively. 

The results for rapidity spectra of secondaries are more model dependent. 
It is possible that the effects of inclusive density saturations are observed at 
CERN SppS energies, and the new RHIC data also give evidences for such 
effects. 
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Figure captions 

Fig. 1. The optical approximation diagrams for the scattering of two 
nuclei. 

Fig. 2. The simplest corrections to the optical approximation which 
correspond (together with the diagrams of Fig. 1) to the rigid target approx- 
imation for the scattering of two nuclei. 

Fig. 3. The distribution on the stripping electrical charge of the projectile 
in C—Ta interactions at 4 GeV per nucleon and its description by rigid target 



approximation [|12 



Fig. 4. Predictions for the associated multiplicity distribution for W^ 
and Z^ events (round black points) together with the experimental data 
{W^ cross points and Z^ white stars) and the minimum bias distribution 
(squares). 

Fig. 5. The x^-spectra of secondary vr"*" (a), tt^ (b), K^ (c) and K~ (d), 
produced in pp interactions at 100 and 175 GeV/c [^ and 400 GeV/c [^]. 
The curves are the Quark-Gluon String Model predictions. 

Fig. 6. The rapidity spectra of secondary 7r+, tt" (a), K^ , K~ (b) and p, 
p (c) produced in Ph + Ph interactions at 158 GeV/c per nucleon |^. The 
curves are the Quark-Gluon String Model predictions. 



13 



References 

[1] Yu.M.Shabelski. Yad.Fiz. 51 (1990) 1396. 

D.Krpic and Yu.M.Shabelski. Yad.Fiz. 52 (1990) 766; Z.Phys. C48 
(1990) 483. 

[2] H.L.Bradt and B.Peters. Phys.Rev. 77 (1950) 54. 

[3] Yu.M.Shabelski. Acta Phys. Polonica BIO (1979) 1049. 

[4] V.M.Braun and Yu.M.Shabelski. Int. J. of Mod.Phys. A3 (1988) 2417. 

[5] V.N.Gribov. ZhETP 56 (1968) 892; 57 (1969) 1306. 

[6] O.V.Kancheli and S.G.Matinyan. Yad.Fiz. 11 (1970) 1305. 

[7] W.Czyz and L.G.Maximon. Ann.Phys. 52 (1969) 59. 

[8] V.Franco and G.K.Varma. Phys.Rev. C18 (1978) 349. 

[9] Yu.M.Shabelski. Yad.Fiz. 47 (1988) 1612. 

[10] G.D.Alkhazov et al. Nucl.Phys. A220 (1977) 365. 

[11] R.D.Viollier and E.Turtschi. Ann.Phys. 124 (1980) 290. 

[12] A.P.Gasparyan, A.P.Cheplakov and Yu.M.Shabelski. Yad.Fiz. 34 (1981) 
1328. 

[13] C.Pajares and A.V.Ramallo. Phys.Rev. C16 (1985) 2800. 

[14] I.V.Andreev and A.V.Chernov. Yad.Fiz. 28 (1978) 477. 

[15] I.V.Andreev and L.A.Khein. Yad.Fiz. 28 (1978) 1499. 

[16] A.S.Pak et al. Yad.Fiz. 30 (1979) 102, 343. 

[17] M.A.Braun. Yad.Fiz. 45 (1987) 1525. 

[18] K.G.Boreskov and A.B.Kaidalov. Yad.Fiz. 48 (1988) 575. 

[19] A.M.Zadorozhnyj, V.V.Uzhinsky and S.Yu.Shmakov. Yad.Fiz. 39 (1984) 
1165. 

[20] V.Franco and R.J.Glauber. Phys.Rev. 142 (1966) 1195. 

[21] Yu.M.Shabelski and D.Treleani. Eur.Phys.J. A2 (1998) 275. 



14 



[22] V.A.Abramovsy, V.N.Gribov and O.V.Kancheli. Yad.Fiz. 18 (1973) 595. 

[23] Yu.M.Shabelski. Nucl.Phys. B132 (1978) 491. 

[24] A.Bialas, M.Bleszynski and W.Czyz. Nucl.Phys. Bill (1976) 461. 

[25] C.Pajares and Yu.M.Shabelski. Phys.Atom.Nucl. 63 (2000) 908; f^- 



ph/9811214 



[26] J.Dias de Deus and R.Ugoccioni. Phys.Lett. B494 (2000) 53. 

[27] A.Capella, A.Kaidalov and J. Tran Thanh Van. Heavy Ion Phys. 9 
(1999) 169. 

[28] Yu.M.Shabelski. Z.Phys. C57 (1993) 409. 

[29] M.G.Ryskin and Yu.M.Shabelski. Z.Phys. C48 (1990) 463. 

[30] G.J.Alner et al. Z.Phys. C33 (1986) 1. 

[31] S.Backovich, D.Krpic and Yu.M.Shabelski. Z.Phys. C53 (1992) 613. 

[32] Yu.M.Shabelski. Sov.J.Nucl.Phys. 50 (1989) 149. 

[33] J.Dias de Deus, C.Pajares and C.A.Salgado. Phys.Lett. B407 (1997) 335. 

[34] J.Dias de Deus, C.Pajares and C.A.Salgado. Phys.Lett. B408 (1997) 417. 

[35] J.Dias de Deus, C.Pajares and C.A.Salgado. Phys.Lett. B409 (1997) 474. 

[36] J.Dias de Deus and C.Pajares. Phys.Lett. B442 (1998) 395. 

[37] A.B.Kaidalov and K.A.Ter-Martirosyan. Yad.Fiz. 39 (1984) 1545; 40 
(1984) 211. 

[38] A.B.Kaidalov and O.I.Piskunova. Yad.Fiz. 43 (1986) 1545. 

[39] Yu.M.Shabelski. Yad.Fiz. 44 (1986) 186. 

[40] A.B.Kaidalov, K.A.Ter-Martirosyan and Yu.M.Shabelski. Yad.Fiz. 43 
(1986) 1282. 

[41] Yu.M.Shabelski. Z.Phys. C38 (1988) 569. 

[42] A.E.Brenner et al. Phys.Rev. D26 (1982) 1497. 

[43] M.Aguilar-Benites et al., LEBC-EHS Coll. Z.Phys. C50 (1991) 405. 



15 



[44] G.H.Arakelyan et al. |hep-pli/010333l . 



[45] F.Sikler et al., NA49 Coll. Nucl.Pliys. A661 (1999) 45c. 
[46] A.Capella and C.A.Salgado. Phys.Rev. C60 (1999) 054906. 
[47] A.Capella and A.Kaidalov. Heavy Ion Physics 9 (1999) 169. 
[48] N.Armesto, C.Pajares and D.Sousa. |hep-pli/0104269| . 



16 



Acp cpQ cpocp poocp 



B i 
a 



6 6 

b 



o o o 
c 



6 6 6 6 

d 




Q 



6 





g 



0.4 



N 



0.35 



0.3 



0.25 



0.2 



0.15 



0. 




0.05 







3 4 

Fig. 3 



Zst 



□ 



n 



n 



^ , « B- 

- □ ^ 

-2 

^ 



-3 
r- 



D ^— ^^; 



n 



i « 



n • I 
n — ^< 

n 

□ i 

n 

n 

n 



20 40 



rt 



f 



° D ^^ 



)( 



t 



^< 4 >< 



i 



60 80 
Fiq. 4 



100 120 140 



n 



PbPb^ K X , A*158 GeV/c 



5 




1 1,5 



,5 3 



r III Q ^ 




PbPb =^ p,p X , A*158 GeV/c 



P A 



5 



:,_„..•-•--+, 



+ 



p T one 



:t- 



X 



±r± 



1 1,5 2 

Fia. 6c 



3 3 



,2( 



^ 7f X , A*1 



120 



2C 



n ▲ 




Fia 



E 




X3 



X 

T3 



b 

T3 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 5c 



E 



b 



400 Ge\ 



-1 




,3 



5,4 



17 



X3 



X 

T3 



b 

T3 



100 and 175 GeV/c 




Fig. 5a 



